Fluctuation AC conductivity of disordered superconducting films 
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We study temperature T and frequency ui dependence of the in-plane fluctuation conductivity 
of a disordered superconducting film above the critical temperature. Our calculation is based on 
the nonlinear sigma model within the Keldysh technique. The fluctuation contributions of different 
physical origin are found and analyzed in a wide frequency range. In the low-frequency range, to <T, 
we reproduce the known leading terms and find additional subleading ones in the Aslamazov-Larkin 
and the Maki-Thompson contributions to the ac conductivity. We also calculate the density of 
states ac correction. In the dc case these contributions logarithmically depend on the Ginzburg- 
Landau rate and are considerably smaller that the leading ones. However, in the ac case an external 
finite-frequency electromagnetic field strongly suppresses the known Alsamazov-Larkin and Maki- 
Thompson ac contributions, while the corresponding new terms and the density of states contribution 
are weakly suppressed and therefore become of relevance at finite frequencies. Moreover, we study 
the conductivity at intermediate and high frequencies, uj > T. 

PACS numbers: 74.40.-n, 74.25.F- 



I. INTRODUCTION 



The behavior of conductivity of a superconductor as 
function of temperature in the vicinity of the supercon- 
ducting transition is governed by superconducting fluc- 
tuations and implicitly contains microscopic and macro- 
scopic parameters characterizing the material. This mo- 
tivated the invariable importance of fluctuations as one 
of the main tools of superconducting studiesEl. One can 
identify three major processes contributing to fluctuation 
conductivity. Fluctuation-induced Cooper pairs, having 
the finite lifetime in the normal domains give rise to the 
increase of conductivity; this is the so-called Aslamazov- 
Larkin (AL) contributioncl. More involved contribution 
arises from the Andreev reflection of quasiparticles from 
superconducting fluctuations accounting for the inter- 
ference due to quasi-particle scattering on impurities; 
this is so-called Maki-Thompson (MT) contributionEHl 
The third fluctuation contribution to conductivity, the 
negative one, is due to the suppression of the quasi- 
particle density of states (DOS) at the Fermi levelEl. 
Impact of all these processes on the dc conductivity 
has been investigated in detail close to the transition 
and far-from itlil. Recently, the extensive studies of 
superconductor-insulator transitional and quest for the 
explanation of mechanisms of high temperature super- 
conductivity-Jjiave revitalized the interest to physics of 
fluctuationsQu, and enriched our understanding of the dc 



fluctuation conductivity!]'!! In this work motivated by 
recent progress in experimental studies of the ac proper- 
ties of disordered superconductorsEJliJ, we address both 
temperature and frequency dependence of the fluctua- 
tion in-plane ac conductivity of films in a wide frequency 
range. The calculations are carried out in the framework 
of the nonlinear sigma model for disordered|-s,unercon- 
ductors via employing the Keldysh techniquetStLS. This 



approach describes the low-energy physics, and is valid at 
energy scales smaller that the quasiparticle elastic scat- 
tering rate and appeal-adequate for the analysis of the 
fluctuation transportOc3. 

The ac in-plane conductivity of superconducting films 
was studied in Refs. |2l],^2] using Matsubara diagrammatic 
technique and the low-frequency, uj <C T, in-plane AL 
and MT contributions were obtained in the leading order. 
The approach chosen in our work has an advantage of be- 
ing more physically transparent than the more traditional 
Matsubara technique since it is formulated in the real- 
time representation and as such does not involve the an- 
alytic continuation. Moreover it applies naturally to out- 
of-equilibrium situationsliMcl Using Keldysh technique 
we obtain the ac conductivity in a wide frequency range. 
We re-derive the well- know n leading low-frequency AL 
and MT contributionsOt3 which are characterized by 
the power law dependence on the Ginzburg-Landau re- 
laxation time in the dc limit and obtain the additional, 
subleading AL and MT contributions to the ac conduc- 
tivity. In the dc limit these subleading corrections behave 
as the DOS contribution, i.e., logarithmically depend on 
the Ginzburg-Landau time. These subleading terms, be- 
ing of little practical importance in the dc limit, become 
of relevance in the ac case, since the leading AL and 
MT terms are strongly suppressed by the applied finite- 
frequency field, while the subleading ones are only weakly 
suppressed. We also calculate the DOS contribution to 
the ac conductivity and find the similar scenario. Al- 
though in the dc case the DOS contribution is sublead- 
ing, it depends weakly on the external-field frequency ui 
and at intermediate frequencies uj ~ T it becomes of a 
similar order as the AL and the MT contributions. 

Recent studies that have been addressing the behav- 
ior of dc fluctuation conductivity appear to diverge in 
details of the final results (in the subleading terms). 
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From this viewpoint an alternative approach based on 
the Keldysh technique is a step towards the understand- 
ing of the possible origin of the existing discrepances. In 
the zero magnetic field, the subleading terms add up to 
C£(3)e 2 In (Tt gl )/ (ir 4 df), where df is the film thickness, 
Tql denotes the Ginzburg-Landau time and C is a nu- 
merical coefficient. The recent paperta finds the result 
C = —231/8 using Matsubara diagrammatic technique. 
A different result with C = —21, is reported in another 
recent studyu. There, the calculation was based on the 
Usadel equation and it is very closely related to the non- 
linear sigma model approach. The results of Refs. || and ^ 
for the individual subleading fluctuation contributions 
agree only on the DOS correction. As a special case of 
our results, at zero frequency we obtain C = —21/2 and 
an additional term depending logarithmically on the ratio 
between the phase coherence rate and the temperature. 
The comparison and discussion of the expressions for the 
individual fluctuation contributions is given in detail in 



Sec. VI 



Apart from the two mentioned papers, there are many 
other older studies of the dc conductivity, based on the 
Matsubara technique, which do not include all the sub- 
leading terms. We do not plan on involving ourselves into 
a comprehensive review of all the interesting works and 
the detailed comparison of all the subtleties since the 
classification of diagrams yielding fluctuation contribu- 
tions to the conductivity in Matsubara technique differs 
from the one within the approach based on the nonlinear 
sigma model. Thus only when all the contributions are 
summed up, the results must coincide and can be com- 
pared. We discuss these issues in more details below, see 



Sec. [VI. 

The paper is organized as follows. In Sec. [Fj 
troduce the model and basic equations. In Sec. 



we m- 



III, the 



MT fluctuation contribution to the in-plane conductivity 
is analysed. In Sees. |^ and ^ we consider the AL and 
DOS contributions, respectively. In Sec. VI, we discuss 



different methods and compare our findings with the re- 
sults k nown in literature. The summary of results is given 
in Sec. VII. In Appendices |X| [FJ, and |c] we present ex- 
pressions for MT, AL and DOS contributions to the con- 
ductivity for arbitrary frequency and Ginzburg-Landau 
time, respectively. 



II. MODEL AND BASIC EQUATIONS 

In this section we introduce the model and basic equa- 
tions that will be used in the rest of the paper. We 
consider a disordered superconductor described by the 
Hamiltonian H = Ha + H mt . The single-particle Hamil- 
tonian in the coherent state basis reads as 



H = / dr i\) a 



(V - ieA) 
2m 



+ Udi S + e<f> 



4a- (1) 



Hereafter h = c = fee = 1- The fields A, <j> and £/di s are 
the vector, scalar and disorder potentials, respectively. 



The electron charge is denoted by e. The spin variable 
is denoted by a =t,4, and the summation over the spin 
indices is implicitly assumed. The interaction is given by 
the Bardeen-Cooper-Schrieffer (BCS) Hamiltonian 



Him = -~ / dr ip^tp^, 



(2) 



where the coupling constant A is positive. The disorder 
potential originates from quenched impurities and it is 
short-ranged. We assume that it is Gaussian distributed 
with the correlator 
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(?7 dis (r)[/ dis (r')} = tf(r-r') 



(3) 



Here v is the bare single particle density of states at the 
Fermi level per one spin projection and r is the elastic 
scattering time. 

In the following we use the Keldysh technique, that 
allows us to perform the disorder average without in- 
troducing replica fields. Then we carry out the standard 
decoupling of the four-fermion terms in the action via the 
Hubbard-Stratonovich fields Q and A, and after integrat- 
ing out the degrees of freedom with the energies higher 
than the elastic scattering rate, we obtain the Keldysh 
nonlineap-,sigma modeLca The partition function takes 
the forrrMtJ 



Z= JvQ 2?Aexp{z5[Q K ,A K ]}, 



(4) 



where the nonlinear sigma model action S is the effective 
action describing the low-energy physics at energy scales 
smaller than the elastic scattering rate r _1 . It holds 
in the dirty limit where the motion of electrons forming 
fluctuation Cooper pairs is diffusive, i.e., the lifetime of 
Cooper pairs is much greater than the elastic scattering 
time. The nonlinear sigma model action reads as 



S[Qk, A/c] = 5a + + Sq, 
where the different contributions are 



(5) 



5a 



v 
2A 



Tr[A^rA K ], 



e 2 v 



Trfcy^jc], (6) 



S Q =—Ti[D(d r Q K ) 2 - AmQK - UefcQ K 
+ 4iA K Q K ]. 



(7) 



The matrix field Q satisfies Q 2 = 1. Here D = TVp/d 
is the diffusion coefficient that carries information about 
the disorder, vp is the Fermi velocity, and d denotes the 
system dimensionality. The check-mark v above the field 
variables indicates that they arc defined in the space that 
is the tensor product of the Keldysh and Nambu spaces. 
The former and the latter arc spanned by the Pauli matri- 
ces &i and Tj, i G {0, x, y, z}, respectively, and we define 
Y = & x (g) fo and 5 = &q ® f z . One uses different nota- 
tion for the same matrices eh = fi for convenience, and 
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do = diag(l, 1). Multiplication in the time-space is im- 
plicitly assumed, and "Tr" includes an integration over 
real space. The subscript K, denotes the gauge trans- 
formed fields: 



K = (k cl 



4>k = 9 ~ &tfc, 
k K = A + V/C, 

do + k q a x ) <g> tq- 



(8) 
(9) 
(10) 



The fields A and <fi are defined in same way as K., while 

A= (A% + A 9 ^) <g> f+ - H.c., (11) 
Ajc(r,t) = e *"^('.*)Ae -<eS *( p '*i. (12) 

The field Qk, is defined in the same way. We have also 
defined f± = (f x ±ifj,)/2. The quantum (q) and classical 
(cl) components of the fields are defined as the half-sum 
and the half-difference of the field values at the lower 
and the upper branches of the Keldysh time-contour, re- 
spectively. The field A cl becomes the superconducting 
order parameter at the mean-field (saddle-point) level, 
while the saddle point equation for Q produces the Us- 
adel quasiclassical equations, where Q plays the role of 
the quasiclassical Greens function. The covariant spatial 
derivative is given by 



d r Qic = V r <5/c - ie[3A/c, Q K ]. 



(13) 



In the absence of the BCS interaction (A = 0), the 



metallic^ 
reads asEi 



point equation for Q obtained from Eq. (^) 



A = UAqW 



A = a z ® t z 
/ 5(t -t'-O) 



(14) 



Ft A*) 



(15) 



After Wigner transforming F t j>(r) we obtain F e (r,t) 
which can be related to the quasiparticle distribution 



function / e (r, t) as F e (r,t) = 1 — 2/ e (r, t). The strat- 
egy is to consider the massless fluctuations around the 
metallic saddle point solution. They can be conveniently 
parameterized aslIS 



Q K {v) = e- W ^' 2 K{v) e 



W{r)/2 



W = UWU~\ (16) 



W = 



WT + — W*T_ Wo TO + W Z T Z 
WqTq + W Z T Z WT + — W*T_ 



(17) 



such that WA + AW = 0. Here we introduced four 
real fields w" t ,(r),w" t ,(r) with a = 0, z representing dif- 
fuson degrees of freedom and the two complex fields 
Wu> (r) , Wtv ( r ) for Cooperon degrees of freedom. The as- 
terisk denotes the complex conjugation. 

In the following, we calculate in-plane ac conductivity 
of a thin disordered film in the normal state but in the 
close vicinity of the transition to the superconducting 
state, T > T c . It can be found using the relation 



Ai=0,A cI =0 



(18) 



where Z is the Keldysh partition function, Eq. (Q). Here 
the index x denotes the x-component of the vector field. 
Since we are neither interested in the weak localiza- 
tion correction nor in the Altshuler-Aronov type correc- 
tions, but in contributions to the ac conductivity caused 
by fluctuations of the superconducting order parameter 
A(r, t), in the following we consider only Cooperon de- 
grees of freedom. We are interested in massless fluctua- 
tions around the metallic saddle point Eq. (|l4]) up to the 
second order in VV, Eq. (|l7j). This means that the film is 
not too close to the transition such that T — T c > G{F C 
where G; <C 1 is the Ginzburg number Gi = {vDdf) -1 . 
Here df denotes the film thickness. We consider linear 
response below, and then we can use the equilibrium dis- 
tribution function F e (r, t)_~_.tanh (e/(2T)) and calculate 
the correlation functionsE2l,E± 



2i 



({w eue2 {q)w*_ e3 _ ei (q)}) = —S ei 

2i 

((w ei , £2 (q)w* (q))) = — <5 ei 



2i 

(K 1 ,6 3 (q)^_ £3) _ e4 (q))) = —S ei - 



~ L K,1-2 L A,1 




■f F ei i A ,l-2 


[Dq2- 


t(ei + e 2 )] [£»q 2 - *(e 3 


+ £4)] 


~ L K,l-2 L A,l 


-2^R.\-2 — F C2 La,1-2 


— FeJ^R, 1-2 


[r»q 2 + 


i(ei + e 2 )] [Dq 2 +i(e 3 


+ **)] 




Afl,l-2 + Fe 2 ^A,l-2 — 





[£>q 2 + i(ei + e 2 )\ [L>q 2 - z(e 3 + e 4 )] 



(19) 
(20) 
(21) 



Here the average ((. . .)} is with respect to the action 
S given by Eq. (||) and it includes averaging over the 
fluctuations of Q, A cl and A 9 . Also, LRiA,i-j = 



f-L R ^(q, e, — 6j)\ denotes retarded/advanced fluctu- 
ation propagators and L~^ { _^ = L^}(q,ei — ej) is the 
Keldysh propagator. The fluctuation propagators read 
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1 Dq 2 — iui 



AttT 



8Tt gl ' 
(22) 

L^{<1, £*>) = [Z^ 1 (q, _)] * = Z^ 1 (q, -_) , (23) 
^A' 1 (q, w) = coth (^) [L^ 1 (q, _) - L^ 1 (q, _)] , (24) 

where the Ginzburg-Landau rate is defined as 



7T V T 



(25) 



Here we expressed the critical temperature as T c = 
2w£)e 7_1 / A /7r, where wd is the Debye frequency and 7 
is the Euler constant. 



•C 1. At higher frequen- 
Dd the Drude conductivity 



This result is valid for lot 
cies, lot 3> 1, the real pa 
is expected to vanish as___ Re[c_>(a;)] oc lj~ 2 , while 
the imaginary part behaves as Iiii[<7d(w)] oc The 
real part describes the dissipation in the system, i.e., the 
attenuation of the external electromagnetic field in the 
sample, while the imaginary part of conductivity gives 
information about its phase shift. If we know either real 
or imaginary part over a wide frequency range, the other 
one is determined by the Kramers-Kronig relation. 



III. MAKI-THOMPSON AC CONDUCTIVITY 

In this section we start the calculation of the ac con- 
ductivity using Eq. (filf). We remind the reader that all 
our results are valid for frequency Li, and temperature T 
much smaller than the elastic scattering rate t _1 , since 
this is the range of applicability of the nonlinear sigma 
model. Therefore in the absence of the BCS interaction 
(A = 0), we obtain the Drude conductivity to be fre- 
quency independent: 

irvDe 2 /j - - - « 

a D {u>) = -— Tr (A^jA^+e + a x A e A- u+e ) 

= J de (F £+L0 -F t _ u ) 

= 2vDe 2 . (26) 



Next we include the BCS interaction and find the su- 
perconducting fluctuation corrections to the ac conduc- 
tivity (|2^) . In the introduction we already explained dif- 
ferent processes that underlie different fluctuation con- 
tributions. We start with the Maki- Thompson contri- 
bution. Considering massless fluctuations around the 
metallic saddle point Eq. ( |l4| ) up to the second order 
in Cooperon degrees of freedom, and collecting all the 
terms of the type {{ww*}} and ((ww*)) in Eq. (|l8|), we 
find MT contribution: 



<7 MT (w) 



vDe 1 I f 

_____ j de 2 de 3 dq [F(e 2 + u) - F(e 2 )} (q>_ e ,__,_ ea __(q) + ™- 



,(q)™ 



(q)»- 



(27) 



Here and in the following, J dq denotes the summation over the discrete component of the wave vector q normal to 
the film (q z ) and integration over the continuous in-plane components (q x ,Qy)- The former gives the factor 1/df. 

In the following we are interested in a system close to the transition to the superconducting state, such that 
t~£ <C T. Then, we find that the leading term reads as 



le 2 D 1 [ + °° [ + ™ 
CTmtM = - - _ 9 - — I de 2 I de 



An 2 d 



f ' 



qdq 



L K \q, e)|_ A (q, e)\ 2 [F(e 2 + _) - F(e 2 )] 
(Dq 2 + i(2e 2 - e))(Dq 2 - i(2e 2 - e + 2w)) : 



(28) 



where ui > 0. The main contribution in Eq. (28) comes simplicity of presentation, in the following we consider 



from small momenta and frequency: Dq , e ~ r GL <C T. two different limiting cases: i) r GL , uj <C T and ii) T,w> 



Therefore, the strategy is to carry out the expansion in 



In order to obtain an expression valid for arbitrary 



this limit and obtain the scries in powers of Tr GL . For frequency uj, one has to consider also other contributions 
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from Eq. (p7j), not contained in Eq. (psj). We provide 
a complete expression valid at arbitrary t gl and w in 
Appendix |A|. 

We start with the parameter region i). Expanding the 
integrand in Eq. (28) for Dq 2 ,uj, eCTwe get the leading 
order terms 



e 7rx — 2 In (x) 

rr GL - 



2nd j 
TC(3)e 2 

TT 4 df 



1+X 2 

In (Tr GL ) 



(29) 



ttx 3 + x 2 - (3x 2 + l) In (x) + 1 
2(x 2 + l) 2 



. (30) 



where x — wt gl . Indexes (1) and (2) correspond to the 
first and second largest contribution. The first term in 
Eq. ([30]) is obtained with the logarithmic accuracy. Note 
that in Eqs. ( p9|) and i'~ ' 



0|) we assumed t g ^ , uj > £>qf nin 



where r , 



denotes the phase breaking ratetl. In the 
, , , in Eqs. (p9|) and (|o|) would ap- 



1 « rr 1 



case u) <C t 

pear r^" 1 instead of w. We point out that result given 
by Eq. ( p9| ) as well as first term in Eq. (ppf), are ob- 
tained using the lowest order expansion of the fluctua- 
tion propagator for e,Dq 2 <C T: L^- (q, e) = i7r/2 and 
L~^ A (q, e) = — -^(Dq 2 + =p ie). In order to obtain 

the second summand in Eq. ( |30| ) one has to use the next 
order expansion of the fluctuation propagators. 

To summarize, the MT conductivity in the limit ui 
but for arbitrary ratio of uj and , is the sum 



Re[cr A , 



-{x)\ =Re[<7MT 



(*)] 



(*)], 



wr GL 



The leading contribution is given by Re[a 



(1)1 



(31) 



This 



is in agreement with Ref. 22 where the MT con- 
tribution was studied in the region i) and the first 

leading term Refcr^r^)] was found. In the dc case 
when t~1 3> tV 1 we get the following dependence 



on the Ginzburg-Landau time from Eqs. (|29|) and 
©: Re[a MT (w = 0)] = e 2 Tr GL In fo/TaJ/frdf) - 
7C(3)e 2 [In (Tr GL ) + In (t*/t gl )/2] /(n 4 d f ). We stress 
that the term proportional to In (t^/t GIj ) originates from 
the last term in Eq. ( |30| ) and in order to obtain it, on 
has to go beyond the lowest order term in the expansion 
of the propagators for small frequency and momentum. 
Also, it is convenient to rewrite the result as 



Re[a MT (oj = 0)] =— r Tr Glj In (t0/t gl ) 



7C(3)e 2 

2TT i df 



[In (Tt gl ) + c In (Tt*)] 



(32) 



with a numerical coefficient c. From Eq. ( |30| ) follows 
c = 1. However, the previous analysis is focused on rel- 
evant dependence on r GL that can be obtained using the 
expansion, and not on the dependence on r^. Therefore 
we additionally analyzed complete MT expression (|2^), 
both analytically and numerically. The complete analyt- 
ical result is giv en in Appendix |A|. In the dc case it can be 
written as (A8). We find that remaining contributions, 
not taken into account in Eqs. ( p9| ) and ( |30| ) give con- 
tributions of the type In (Tr^) and change the numerical 
coefficient in front of this term from c=ltoc~3.5. As 
expected, they do not influence singular dependence on 

Tgl- 

Now we consider the parameter region uj,T t gl . 
The ratio of temperature T and frequency uj is arbitrary. 
Then, in the similar way as in the case i), we obtain in 
the leading order two contributions 



4 2 / T N 
Re[a$(w)] =- rr - ) In (Tr GL ) 



7T df \U) 



x Re 



1 



Ra[^M]=^(|) tanh(^) 



, (33) 
(34) 



with the logarithmic accuracy. At intermediate and high 



frequencies uj ^ T, the terms Cmt an d of similar 

order with the weak domination of CTmt as the frequency 
increases or the GL rate decreases, while for lower fre- 
quencies uj <C T the leading contribution is given by 
Eq. ([m]). The regions i) and ii) overleap at frequen- 
cies satisfying <C uj <C T '. Notice that there the 
subleading low-frequency term ( p0| ) matches the contri- 
bution (|33j), and also leading low- frequency contribution 
( p9| ) matches the term To conclude, analyzing re- 

gions i) and ii) we see that the leading term in the dc 
and low-frequency range, becomes more suppressed by 
the finite-frequency external electromagnetic field than 
the subleading term (in the low- frequency range), and 
therefore they become of similar magnitude at uj > T. 



r (2) 



IV. ASLAMAZOV-LARKIN AC 
CONDUCTIVITY 



In this section we calculate the Aslamazov-Larkin ac 
conductivity. It is determined by the terms of Eq. ( |is| ) 
which contain four fields that describe Cooperon degrees 
of freedom. We obtain 
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c A lM = - 2(2tt) 8 Zi J deide2de4de5dqidq3glzg3z 

U (F(ei) [tD*(-qi, -e 2 , -ei + w)u)(-qi, ei, e 2 ) - w(qi, e 2 , ei - w)u)*(qi, -ei, -e 2 )] 

+ F(ei - w) [w(qi, e 2 , ei - w)uj*(qi, -ei, -e 2 ) - w*(-qi, -e 2 , -e x + w)ti;(-qi, e lf e 2 )] } 

[w(q 3 ,£5,e4 + w)u;*(q 3 , -e 4 , -e 5 ) - w*(-q 3 , -e 5 , -e 4 - uj)w(-q 3 , e 4 , e 5 ) 

+ w(q 3 ,e5,£4 + w)w*(q 3 , -e4, -e 5 ) - w*(-q 3 , -e 5 , -e 4 - w)w(-q 3 , £4, e 5)])) (35) 

After performing Wick contractions, the summands of the previous expression can be classified with respect to the 
number of the distribution functions they contain: one, two, and three. It turns out that the term that contains 
one distribution function F, exactly nullifies after the energy integration. The reason is that all the zeros of the 
corresponding denominators are in one half-plane of the complex plane. The term with two F gives the leading 
contribution of Eq. ( |35| ) in the vicinity of the transition into superconducting state. After exact integrations it reads 
as 

°® =T^T / + °° ^ r ^ wwr?^- M [G(z)][G{z) - G(z - 2ia)] 
2ir 2 a 6 d f J_ 00 J Q \G(z)\ 2 G(z - ia) 

[G(z*) + G(z - ia) - G(z* - ia) - G(z - 2m)]. (36) 



Here we have introduced the short-hand notation 
G(z) = ^(1/2 + z)- ^(1/2) + '' 



4nT' 

z = x + iy, z* = x — iy. 



8?W ^ 
(38) 
(39) 



Next we evaluate expression (|36|) analytically for small 
frequencies, while it can be evaluated numerically at ar- 
bitrary frequency. In the limit , w <C T but for arbi- 
trary ratio of co and 



•GL! We nIld 



Re [a 



(1,1)1 

AL 



Reki 1 ^] 



26 (T ^ 1 
" j v t glJ 2 2 

-HI+^tIJA)] 
7C(3)e 2 1 



V M 'GL 

,2_2 



[wt gi arctan(ojr GL /2) 

GL 

(40) 

4wt G l arctan(LLir GL /2) 



+ 21x1(1 + ^/4) 

1 + ^/4 



2 2 i 
uj In • 



T 2 r 2 

'GL 



(41) 



The index (1,1) denotes the leading term and (1,2) the 
subleading term. Notice that in order to obtain the result 
( fi"T| ) one has to go beyond the lowest order expansion 
of the fluctuation propagators for small frequency and 
momenta. The last term in Eq. ^l] is determined with 
the logarithmic accuracy. 

Next we consider the term where the distribution func- 
tion appears three times. Its full dependence for arbi- 
trary r GL and ui is given in Appendix Here we calcu- 
late its leading contribution for , w -C T. It is of the 



same order as cr A L 2 \ and reads as 



!^-lRe[J 2 ]Im[7i] 



7T at lu Jo 

x {2Rc[L A (q, e)L A (q, e ~ 
+ L R (q,e)L A (q,e~uj)}. 



dqq 3 



+00 



de 



(42) 



In this equation L A m denote the lowest order in the ex- 
pansion of the fluctuation propagators for small frequen- 
cies and momenta. The functions I\ and J 2 are defined 
as 



h=- 
h=- 



tanh(e/2T) _ . tt 
(2e-z0) 2 ~ ~*8T' 

7C(3) 



+00 dg tanh 2 (e/2T) 



(2e-i0) 2 2tt 2 T' 
After evaluation of expression (fi"2|), we obtain 

7C(3)e 2 1 



(43) 
(44) 



2vr 4 d 



/ wr GL 



■ 4 arctan(wr GL /2) 



In 



I+^gl/4 
T 2 t 2 



(45) 



with the logarithmic accuracy. 

To summarize, in the limit t^,lu <C T we find the 
leading term given by Eq. (|i(]) and the next leading 

order contribution is the sum of ctXl, 2 ^ and a 1 ^ 1 ^ given by 
Eqs. (Ell) and respectively. They together give 
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Re[cr AL (o;)] = 



2e 2 T 



arctan 



wdf ui 
28C(3)e 2 



wt gl 



■In 1 



arctan 



■In 1 



■In 



l + ^ L /4 



(46) 



The first leading term is in the agreement with calcula- 
tion of Refs. |2l| , |22| , p7| , while the second one is new. From 
Eq. ([46]), in the dc case we obtain 



Re[a AL (0)] = — -Tt, 



2nd 



f 



14C(3)e 
n*d f 



■ In (Tr GL ). (47) 



Apart from the well known power law singularity of the 
AL dc contribution, we get the additional logarithmic 
divergence. Moreover, from Eq. (^6|), one gets that for 
uj 3> Tq L , the subleading contribution is logarithmically 
suppressed with increasing uj, while the leading term de- 
cays much faster as w" 1 . So, although in the dc case 



the subleading term <t A l' 2 ' ) + a < ^, 1 " > could be considered as 
negligibly small, in the ac case their sum becomes more 
relevant. 



V. DENSITY OF STATES AC CONDUCTIVITY 

In this section we calculate the density of states contri- 
bution to the conductivity. We start form Eq. (18) and 
collecting all the terms of the type {{ww*)) and ((ww*)) 
we find 



Cdos(w) 



vDe 2 If ( 

16^3 uj / de 3 de 4dq|i ? (e 3 + w)((w £3ie4 (q)w* £4 _ e3 (q))) - F(e 3 - w)((w e3i£4 (q)w* £4 _ £3 (q))) 



F(^)((We 3 ,eAci)w* u ^ e3 _Jq))) + F(e 3 )((u} £3 , £4 (q)™* 



(q)))} 



(48) 



In Appendix ^|, we provide an analytic form for arbitrary frequency and arbitrary Ginzburg-Landau time, while here 
we focus on the case when the system is close to the superconducting transition. Then, the leading contribution comes 
from the term 



CTdOsM = + 



iDe 2 1 f +oc 



An 2 df 



de 3 de J AqqL K \^e)\L A {c i ,e)\ 2 [ 
tanh (||t) 



tanh(^) 



tanh (^) 



[Dq 2 + i(2e S - e)] 2 [Dq 2 - i(2e 3 - e)} 2 
tanh(^) 



[Dq 2 + i(2e 3 - e)}[Dq 2 + i(2e 3 - e - 2u)\ [Dq 2 - i(2e 3 - e)][Dq 2 - t(2e 3 - e + 2w)] . 



We obtain the leading contribution for t g ^ «C T, uj to the ac conductivity to be 



Re [o- D osM] 



2e 2 T 



In (TV OL )< - Im 



TT 3 df UJ 

From this equation we get the limiting cases: 



2nT) 



T 

27r-Re 

UJ 



' ' 2 % 2nT 



(49) 



(50) 



Re[cr D osM] 



ln(Tr GL ), r- 1 « W «T 
In (Tr GL ) (T/uj) 2 In (uj/T), « T « uj. 



(51) 



r 



Analyzing Eq. (Q) in the limit uj <C -C T, we find This result is also valid in the dc case, and we see that 

the DOS contribution is almost frequency independent at 
21e 2 £(3) uj C T. Characteristic frequency for the DOS contribu- 

<n±df ^ ^ Tgl '' tion is determined by temperature T. This is in contrast 



Re[a DOS (w)] 
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with the AL and the MT corrections whose character- 
istic frequency is determined by T — T c . They decrease 
as power law ~ w _1 with frequency, for r^ 1 <u<T. 
Therefore, at cj ~ T the DOS becomes of the same mag- 
nitude as the MT and the AL contributions. 

Notice that the subleading contributions in other fluc- 
tuation corrections are of the same order as the lead- 
ing DOS contribution at low frequencies. Also, notice 
that the DOS conductivity logarithmically diverges for 
Tq L — » 0, even for finite frequency u > 0, and it is nega- 
tive. We remind the reader that all our results are valid 
for frequency, and temperature much smaller than the 
elastic scattering rate, since this is the range of applica- 
bility of the nonlinear sigma model. 



VI. DISCUSSIONS 



We are now equipped to compare different methods 
used for the analysis of the fluctuation transport and dis- 
cuss some of the results present in the literature. First 
we would like to point out that classification of the dia- 
grams into the different fluctuation contributions to con- 
ductivity in the conventional Matsubara diagrammatic 
technique^ differs from that in the approach based on the 
nonlinear sigma model, which we employ in our work. In 
the latter, the physical origin of different contributions 
appears to be more transparent, as we discuss below. 
However, when comparing our and the conventional dia- 
grammatic approaches, the final results summing up all 
the contributions in both calculations have to coincide. 
More specifically, we expect that sum of our MT and 
DOS corrections should be equal to the sum of the MT 
and DOS corrections in the diagrammatic approach, im- 
plying that also the AL correction is the same in both ap- 
proaches. The AL correction contains two contractions 
of the Cooperon fields and they correspond to two fluctu- 
ation propagators appearing in the diagrams for the AL 
correction in the conventional approach. 

Having said that, we turn to our classification of the 
fluctuation contributions. All the terms in Eq. ([l8]) that 
contain a convolution of w and w* as well as w* and w 
give the MT contribution, while all the terms contain- 
ing a convolution of w and w* as well as of w and w* 
are collected into the DOS contribution to the conduc- 
tivity. After redistribution of terms into the MT and 
the DOS contributions, the remaining ones belong to the 
Alsamazov-Larkin contribution, Eq. (35). These remain- 
ing terms contain four fields and applying the Wick the- 
orem we get products of two contractions. 

Let us analyze in more detail the structure of the terms 
in the DOS contribution. We notice that the first and the 
second term in Eq. have slightly different structure 
than the third and fourth term. Also, the expression for 
the change of the density of states of quasiparticles due 



to superconducting fluctuations reads as 

<M £ ) = -7^3 J de i J dq((«)e > e 1 (q)w!L £1 _ e (q) 

+^ ei (q) W * £1: _ e (q)», (53) 

and then the first and the second term in Eq. ((Iq) in the 
dc case can be rewritten as 2De 2 J deSi/(e)d e F(e), while 
the third and fourth term have different structure. How- 
ever, since all the contributions include the same type 
of the Cooperon degrees of freedom, we group them into 
the DOS correction. 

Note that the important progress in treating fluctu- 
ations within the Keldysh nonlinear sigma model was 
achieved in Refs. [l7 28, where the leading fluctuation 
contributions to the dc conductivity where calculated. 
However, there the contributions to the conductivity cor- 
responding to the two last terms in Eq. (^) were missed. 
They obtained the factor —7 for the DOS correction. As 
follows from our consideration, in the dc case the first two 
summands in Eq. ( [f8| ) give numerical factor —14 and the 
last two give —7, and they together give factor —21 in 
Eq. (^2|) [note that there had been a misprint in the pub- 
lication of the present authors, Ref. |2^, in Eq. (58), where 
the numerical factor —21 should stand instead of —7]. 

Next we summarize the results on the ac in-plane 
transport known in the literature and compare them with 
our findings. In Ref. ^l] the leading ac low-frequency AL 
term was found and later in Ref. 22 the leading low- 
frequency MT contribution was calculated. These re- 
sults are in an agreement with our results. In Ref. |27] 
the ac DOS correction was obtained. All these stud- 
ies employed Matsubara diagrammatic technique. As we 
explained above, the sum of the DOS and the MT con- 
tributions obtained in this approach should be compared 
with our results. However, because only the leading MT 
contribution was calculated in Ref. ^ (it corresponds to 
the leading part of the so-called anomalous contribution) 
and the subleading MT terms that are of the similar mag- 
nitude as the DOS contribution were not considered, we 
can not compare our result for the DOS correction with 
the findings of Ref. p7j 

We now analyze the dc case where the subleading cor- 
rections in all the fluctuation contributions were consid- 
ered, contrary to the ac case, allowing us to undertake a 
more detailed comparison. We have found in the previous 
sections that close to the superconducting transition 



Re[cr(w = 0)] 



2nd f 
, e 2 C(3) 



Tr GL [1 + 2 In (t /t g 



In (Tt gl ) (c al + Cut + Cdos) 

(If 7T 



e 2 7C(3) 
d f 2tt 4 



cln (Tt^) 



(54) 



for 3> t gl . Here Cj, where i takes values al,mt, dos 
describes contributions of different physical origin, the 
AL, MT and DOS. We find c AL = 14, c MT = -7/2 and 
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Cdos = —21. The term In (Tt^) comes from the MT 
contribution and c « 3.5. The leading terms arepgiven in 
the first line of Eq. (^3) and are well established!!], unlike 
the remaining subleading terms. 



The recent publication, Ref. ||, studied the dc conduc- 
tivity using both the Usadel equation and the Matsub- 
ara technique. The results based on the Usadel equation 
give the total numerical coefficient C = Cj = —21 in 
front of the logaritmically singular term In (TY GL ). The 
authors of Ref. || used different gauge where the elec- 
tric field is given by E = — V0, and found c MT = —14, 
c D os = —7 and c AL = 0. In order to compare the results 
based on the Usadel equation of Ref. |9| with ours, one has 
to consider the sum of the DOS and the AL terms, since 
the expressions for the MT dc contribution, our Eq. (I^t) 
and Eq. (78) in Ref. | for zero magnetic field, coincideea. 
For the sum of the AL and the DOS terms there is full 
agreement, while the ^efficients c MT are different due to 
inaccurate evaluatiorO of Eq. (78) in Ref. ^. The other 
approach in Ref. ^| based on the conventional diagram- 
matic approach gives the numerical coefficient C = — 21 
as the sum of c DO s = —14, c AL = +14 and c MT = —21. As 
we discussed at the beginning of this section, in order to 
compare our results with the ones based on the conven- 
tional diagrammatic approach, one should compare the 
AL terms and the sum of the MT and the DOS terms. 
The result of Ref. || for the AL contribution agrees with 
our result. However, there is a disagreement when com- 
paring the sum of the subleading DOS and MT contri- 
butions of Ref. 0. Since Eq. (|A8j ) also appears as a part 
of the MT contribution in the diagrammatic approach, 
the above mentioned inaccurate evaluation of this term 
in Ref. || explains the disagreement in the sum of the MT 
and the DOS terms. The two above comparisons of our 
findings with the results of Ref. || provides us with the 
comparison of each individual fluctuation contribution, 
and we infer agreement in the AL and the DOS results, 
while the disagreement in the MT term is resolved. 



Note that in the recent paper of one of the present 
authors, Ref. ||, it was found C = —231/8 using the Mat- 
subara diagrammatic technique. (In Ref. ^ see Table 
1 and region I.) There c AL = —7/8, c MT = —14, and 
Cdos = —14. The fact that the results on the subleading 
AL term of the present work agree with those in Ref. H, 
it adds to the confidence to validity of the latest findings 
that c AL = 14, so the revealing the reasons for the re- 
maining discrepancy in the subleading AL contributions 
in Ref. |^ requires more analysis. It also requires addi- 
tional inspection if there should be an additional sub- 
leading MT contribution in Ref. || that depends logarith- 
mically on the ratio of the phase coherence time and the 
Ginzburg-Landau time that would change the coefficient 



VII. SUMMARY OF RESULTS 

We considered a disordered film above the supercon- 
ducting transition and studied influence of superconduct- 
ing fluctuations on the ac in-plane conductivity. Our ap- 
proach is based on the nonlinear sigma model for dirty 
superconductors that is valid for frequencies smaller than 
the elastic scattering rate. Within this limitation, we 
provide analytical expressions for the ac in-plane conduc- 
tivity for arbitrary frequency u> and arbitrary Ginzburg- 
Landau time r GL . We analyzed them in more detail in 



the following frequency regions: 



,uj <S^T and arbi- 



Cmt — 14 to c^jt 
below Eq. (|3l|). 



-21/2, see discussion in Sec. [II 



trary ratio between ui and r^" L j and ii) T, ui 3> t~ l but 
arbitrary ratio between uj and T. These regions have an 
overlap and provide us with conductivity behavior in the 
wide frequency range. 

We calculated: (a) the Maki- Thompson contributions 
to the ac conductivity given by Eqs. (^9|,[50|) in the re- 
gion i) and by Eqs. (|3^j3^ ) in the region ii); (b) the 
Aslamazov-Larkin contribution, Eq. ( |l6|); (c) the den- 
sity of states contribution given by Eq. (p2|) in the region 
i) and by Eq. ( |50| ) in the region ii). Our results agree 
with the existing leading MT and AL contributions to 
the conductivity, that were known for low frequencies 
w CT, see Eqs. ( p9| ) and (|40|). The obtained expressions 
for the subleading low-frequency behavior of the AL [sec- 
ond line of Eq. (|e|] and the MT [Eq. (|o|)] corrections, 
and intermediate- and high-frequency behavior of the MT 
[Eqs. @ and Q] and the DOS [Eq. |(]] corrections are 
new and provide us with a complete physical picture. For 
example, in the dc case the DOS contribution could be 
ignored with respect to the others. However, at frequen- 
cies uj <C T it is almost frequency-independent in contrast 
to the other two contributions which get strongly sup- 
pressed by the external electromagnetic field. Therefore 
the DOS contribution becomes of the same magnitude as 
compared to the other two ones for ui ~ T . 

A similar approach to the one employed here can be 
used for calculation of the c-axis transport. In the mea- 
surements of the c-axis reflectivity spectra in YBa2Cu40g 
single crystals in Refs. [50 31 it was found that the c-axis 
optical conductivity shows a transition from a Drude-like 
to a pseudo-gap like behavior with the decrease of the 
temperature. In Ref. ^?], the influence of the supercon- 
ducting fluctuations on the c-axis conductivity was con- 
sidered and the non-monotonic frequency behavior that 
corresponds to the pseudo-gap was found. Its origin is 
related to the nonmonotonic frequency behavior of the 
DOS correction. Note that in the case of c-axis transport, 
the positive AL and MT contributions are suppressed by 
the interlayer transmittance. and compete with the neg- 
ative DOS contributioi£3l23. At the intermediate fre- 
quencies, the DOS correction dominates and one obtains 
nonmonotonic behaviorEZI. However, the difference be- 
tween the minimal and the maximal conductivity of this 
nonmonotonic behavior is small and only logaritmically 
increases as the transition is approached. Therefore, it 
would be interesting to calculate subleading AL and MT 
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contributions and see how they influence this result. In 
the c-axis conductivity the leading AL term behaves as 
r,? L and the leading MT term as r GL , while the DOS term 
shows, logarithmic dependence on the Ginzburg-Landau 
timeEj. One can expect that the subleading term in the 
AL contribution might depend linearly on the Ginzburg- 
Landau time and might play a role in the final result as 
the transition is approached. 
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Appendix A: Maki-Thompson contribution 

In this appendix we present the complete expression 
for the MT contribution (|2^), valid for arbitrary t gl and 
arbitrary frequency. After performing one integration, it 
can be written as the sum of the following terms: 



%T 1 



8n z a jo 



coth(27ry) Im[G(z)] 
da; ; , — — {2Re|G(z)J — G(z — 2ia) — G(z — 2ia)\ , 



<y mt2 



-oo ju x-ia \G(z)\ 2 

I r+co poo 

—5- / dy / dx\ + coth(27ry)[2iIm[G(z)] + G(z - 2ia) - G(z* - 2ia)\ 
37Ha J_ oq Jo 

- coth[27r(y + a)][2G(z) - 2Um[G(2i(y + a))} - 2G(z* - 2ia)]} 

Here we have introduced the short-hand notation 

G{z) = yj(l/2 + z)- ^(1/2) + " 



1 

(x — ia)G(z) 



8Tt g 



a 



We used the integral 



dx 



tanh i 



(x + a)(x + b) a — b 



47rT' 

z = x + iy, z* = x — iy. 



*<5-7)-*(H 



and the identity 



- Z) = (-l) n 1p (n) (z) + (-1)™7T 





2z7r tanh a 
a — b 



d n COt 7TZ 



Im(a),Im(6) > 0, 
Im(a) < 0, Im(6) > 



dz n 



(Al) 

(A2) 

(A3) 

(A4) 
(A5) 

(A6) 
(A7) 



for a positive integer number n. Note that er MT i coincides with Eq. (p8|), while the remaining terms in Eq. (|27]) are 
denote d by £t M t2- The main contribution close to the superconducting transition comes from er MT i, as discussed in 
section 



III 



Im 2 [G(z)] 



(A8) 



In the dc case from Eqs. flAl| ) and (A2) we find 

/-f-oo p oo 

-oo dy l ml „ da: ^sinh 2 (27ry)|G(z)| 2 
Note that in the dc case one should regularize these integrals using the phase coherence time t^ 1 which determinesQ 

Xmin ~ T^ 1 /(4irT) . 



Appendix B: AL contribution 

In this section we provide the complete expression for the AL contribution (^) valid for arbitrary frequency and 
arbitrary Ginzburg-Landau time. It can be written as the sum of the two terms, one is given by Eq. (B™ and the 
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other one is given by 

/■+0O 



dy j da;| — coth(27ry — 2tto) 



x [G (z*) -G(z* - ia)} {G (z*) + 2ilm[G(z - ia)} - G(z - 2ia)} 
4ir 2 a s G(z*)G(z - ia) 



+ 



x coth(27ry) {—G {z* - ia) + G(z - ia) - 2Um[G(z)]} 
8TT 2 a 3 \G(z)\ 2 G{z - ia)G (z* - ia) 

x (G (z*) {-G(z - 2ia)G (z* - ia) + 2G(z) [G (z* - ia) + G(z - ia)} - G(z - ia) 2 } 

-G(z) G (z* - ia) 2 + G(z - ia)G (z* - 2ia) ) 

x coth(27ra) 
+ 8ir 2 a 3 \G(z)\ 2 G(z + ia)G (z* + ia) 

x G (z*) ([G(z - ia) ~ G(z + ia)} [G{z - ia) - G(z + ia) + G(z + 2ia)}G (z* + ia) 
- G(z) {G(z + ia) [G (z* - ia) - 2G (z* + ia)} + \G(z - ia)\ 2 } ) 

+ G(z)G(z + ia) [G (z* + ia) - G (z* ~ ia)} [-G (z* - ia) + G (z* + ia) - G (z* + 2ia)} } 



(Bl) 



The main contribution from this term when the system is close to the transition is given in Eq. (|42| 
In the dc case we find that total AL contribution for arbitrary r GL can be written as 



/ + oo rco 
dy / dxl - 
-ao JO 



xlm 2 [G'[z}} 



Att 3 \G{z)\ 
4coth(27n/){ 



|G(z)| 2 sinh 2 (27ry) 
{\G(z)\ 2 Re[G 2 (z*)G'(z)G"'(z)} - 2Re[G 3 (z)[G'(z*)} 2 G"(z*)}} 



Im[G'(z)/G(z*)}\G'(z)\ 2 +Im{[G'(z*)} 3 /G(z*)} + 2Im[G'(z)G"(z)} 



2|G(*)P 



+ Im • 



G'(z) 
G(z) 



- 21m 



G'(z)G"{z) 



G 2 (z) 



}}• 



(B2) 



This expression contains many terms that are unimportant close to the transition, see Sec. [Tv| . 



Appendix C: DOS contribution 



In this section we provide the complete expression for the DOS contribution (48), valid for arbitrary frequency and 
arbitrary Ginzburg-Landau time. It can be written as the sum of the two terms: 



1 



CdosI 



0DOS2 



47r 2 a 2 
i 



dx ^^h miGiz)} [G{z*) ~ G(z* - 2ia) + ia[G'(z - 2ia) + G'(z* - 2ia)}} , (CI) 



87r 2 a 2 



dy / dx 

-oo JO 



\G(z)\ 2 

2ilm[G{z)} + G(z* - 2ia) - G(z - 2ia) 
W) 



+ 2ia 



G'(z - ia) - G'(z* - ia) 



G(z — ia) 



coth(27ra) 

W) 



[2G(z) - G(z - 2ia) - G{z + 2ia) - 2iaG'(z) + 2iaG'(z - 2ia)] 



coth(27ry) 
(C2) 



we used the same notation as in the previous appendix, Eqs. ( A3 ) , ( |A4| ) , and (A.5). Note that here cr DO si coincides 
with Eq. ([l9|), while the remaining terms in Eq. ( ^ ) are denoted by cr DO s2- The main contributio n c lose to the 
superconducting transition comes from <r DOS i, as discussed in section |v|. In the dc case, from Eqs. (31) and ( |C2| ) 
follows 



/+oo poo 
dy / dxl - Re 
-oo Jo 



G"'(z) 



4tt 3 G(z) 



7n 2 \G(z)\ 2 { A MG{z)}lHG"{z)} + 6lm[G(z)}Re[G"(z)} 



2lm[G l (z)]Rc[G 2 {z*)G'(z)} 
\G(z)\ 2 



}}■ 

(C3) 
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However, it turns out that many terms in this general expression are unimportant close to the superconducting 
transition, see Sec. |v|. 
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